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We study, using ąuantum Monte-Carlo simulations, the bosonic Kondo-Hubbard model in a two 
dimensional sąuare lattice. We explore the phase diagram and analyse the mobility of particles 
and magnetic properties. At unit filling, the transition from a paramagnetic Mott insulator to a 
ferromagnetic superfłuid appears continuous, contrary to what was predicted witli mean field. For 
double occupation per site, both the Mott insulating and superfłuid phases are ferromagnetic and 
the transition is still continuous. Multiband tight binding Hamiltonians can be realized in optical 
lattice experiments, which offer not only the possibility of tuning the different energy scales over 
wide ranges, but also the option of loading the system with either fermionie or bosonic atoms. 
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I. INTRODUCTION 

In condensed matter Systems, the interaction between 
mobile particles and fixed magnetic impurities, known 
as Kondo physics, has been a very important topie for 
the last 50 years. From the original explanation of the 
resistance minimum in metals with magnetic impurities 
by Kondo, to the investigation of the properties of heavy 
fermion materials, the interaction between particles and 
localized spins has revealed a variety of interesting 
physical phenomena GL 0- Indeed, the competition 
between magnetic ordering and singlet formation in 
Kondo and related materials has offered some of the most 
fundamental examples of ąuantum phase transitions 0, 
and investigations of the effects of interplay of the distinct 
spin and particie contributions to the susceptibility, and 
of dilution0,0, are at the frontier of the inyestigation of 
many materials, including the ‘115’ heavy fermion family 

m 

In addition to these solid State Systems, with the recent 
experimental advances in ultracold atomie physics, it is 
now possible to build Systems of atoms on optical lattices 
with atoms occupying different bands 01 - This opens 
the possibility to use the atoms located in the lowest 
band as localized particles, magnetic centres, which will 
interact with mobile particles located in higher energy 
bands. Such Systems would be analogues of Kondo 
problems but with the possibility to use bosonic particles 
instead of fermionie ones, Systems that have not been 
extensively studied and are not available in condensed 
matter physics. 

We will study here a system similar to the Kondo- 
Hubbard lattice problem [2|, H0-ll2| with interacting spin 
1/2 bosons instead of fermions. Mobile bosons are free 
to move on the lattice and interact repulsively on site. 
In addition, there is an antiferromagnetic (AF) coupling 
to an ensemble of spin 1/2 magnetic centres, one for 
each site of the lattice. This model was introduced by 
Duan jTfij] to describe the following system: the localized 
bosonic species (‘spins’) are atoms in the lowest band of 


an optical lattice with a potential barrier which prohibits 
tunneling. The mobile species occupy an upper band, 
to which they have been excited through the applicaton 
of periodic Raman pulses, which allows tunneling. This 
model was studied in detali with different analytical 
techniąues by Foss-Feig and Rey 0 - In 0 ], exact 
results were derived for the smali and large interaction 
limits and the intermediate regime was studied using 
mean-field theory, for different densities of particles. For 
one mobile particie per site, they observed a first order 
transition between a Mott insulator (MI) phase and a 
superfłuid (SF) phase as the interaction is lowered. In the 
Mott phase, there are singlets of bosons and spin and no 
long rangę magnetic order whereas the SF phase shows 
long rangę ferromagnetic (FM) order for the bosons and 
the spins. For two or morę particles per site, the Mott 
phase is already ferromagnetic and the transition to the 
ferromagnetic superfłuid State is continuous. 

In this paper, we will use exact ąuantum Monte 
Carlo simulations to study this bosonic Kondo-Hubbard 
model 0 and determine exactly the phase diagram 
and magnetic properties of the system at zero and finite 
temperatures T and compare with results pręyiously 
obtained with mean-field approximations [131. 114l|. In Sec. 
II, we introduce the model, the numerical techniąue we 
used and the ąuantities we will measure to characterize 
the phases. In Sec. III, we study the transport properties 
and Green functions of the system to draw its phase 
diagram at T = 0. In Sec. IV we analyse in morę detail 
the naturę of the ąuantum phase transitions. Sec. V is 
devoted to a careful analysis of the magnetic properties 
in the ground State and Sec. VI presents the evolution 
of the phases observed at T = 0 as the temperaturę is 
inereased. 


II. BOSONIC KONDO-HUBBARD MODEL 

The system we consider includes two types of objects 
that are coupled antiferromagnetically (AF): spin-1/2 
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a: spin 1/2 bosons 
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b: 1/2 spins 


Figurę 1: Schematic representation of the model. Spin-1/2 
bosons (upper layer a) are moving through the lattice with 
hopping i and are subject to on-site repulsion U. They 
interact antiferromagnetically with strength V with the lower 
spin layer b. The lattice represented here is a ID chain; the 
studied system is a 2D sąuare lattice. 


bosons which hop on a 2D sąuare lattice and, on every 
site, a fixed spin-1/2 magnetic impurity. In the following, 
we will denote the bosons as particles of type a and the 
fixed spins as particles of type b. An easy way to visualize 
this system is to use two layers, a and b , one filled with 
bosons, the other with spins (Fig. [[]). The Hamiltonian 
reads 


R = -t J2 ( a Us> + h - c -) 

(1) 

U ^ 

+ ~2 Z__J ( n a,i - 1) 

i 

(2) 

± 4^^ S 0ii • S M 

(3) 
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The operators aj a and cq jCr create or destroy an a-type 
boson of spin er on site i. The system is a 2D sąuare 
lattice with L 2 sites where L is the length of the lattice. 
There is always one b spin per site and we will vary the 
number N a of a bosons. The Hamiltonian includes a 
hopping term (Eq. [T]) and an on site repulsion (Eq. [2]) 

for the a bosons. n ati = na,i,a = a l,a a i,a is tlie 

total number of bosons on site i. The hopping parameter 
t sets the energy scalę and the on-site repulsion energy is 

U. 

The last term (Eq. [3]) of the Hamiltonian is an 
antiferromagnetic coupling between the boson magnetic 
moment and the fixed b spins. S 0) j = i; S* J 

gives the spin of the bosons and its components are given 
by 




(4) 


where the S“ , are the three standard spin 1/2 matrices. 
The b spins are also described by these matrices S b,i = 
(S% SjF, Sfii). To specify a State of the system one 
should give the State S gF = ±1/2 of each b spin and 
the number n a ,i,cr of up or down a bosons present on 
each site i. The conventional ąuantum number s will be 


used when discussing the possible eigenvalues of angular 
momentum S 2 = s(s ± 1). 

The last term of the Hamiltonian is the Kondo 
interaction, here in the form used in the Kondo insulators 
where the moving particles interact with a network of 
magnetic moments. This is different from the original 
Kondo problem where the moving particles are coupled 
to a smali number of magnetic “impurities” distributed 
randomly p} and morę similar to the “Kondo lattice” 
m. Other differences with the original Kondo problem 
are that our moving particles are not free but interacting 
with each other and, of course, they are bosons and not 
fermions. Studying an equivalent spin 1 model would be 
interesting but is morę demanding numerically. The spin 
1/2 model allows us to compare with the results from 
[13Ł Il4| , which also propose an experimental realization 
of the model. Furthermore the qualitative physics should 
not be different with a spin 1 model. 

To study this system, we used the ąuantum Monte 
Carlo SGF algorithm [H, [l7| that allows exact 
calculations of physical observables at finite temperaturę 
on clusters of finite size (up to L = 14). We are especially 
interested in one and two-body Green functions that are 
possible to calculate with the SGF algorithm. To extract 
the properties of the ground State, we used large inverse 
temperatures fi = 1/kT , up to (5t = 25. 

We studied the one-body Green functions for the 
bosons G a ,a{R), 

G a ,*{R) = + a \+R,a a i,c)- ( 5 ) 


The condensed fraction p(k = 0) is the Fourier transform 
at k = 0, p(k = 0) = J2r a G a ,a{R)/L 2 - The superfluid 
density p s can be measured using the standard relation 
with fluctuations of the winding number p s = (W 2 )/4f 0 
as the total number of bosons is conserved [Isj . 

We also studied anticorrelated two-body Green 
functions which describe exchange of particles or spins 
at long distances, which then correspond to opposite, 
anticorrelated movements. They are generally important 
for multispecies Hamiltonians with repulsive interactions 
where exchanges are the dominant effects in the strongly 
interacting regimes [21j. In this case, they are 
conveniently expressed in terms of spin degrees of 
freedom 


G aa {R) — ~b b-C-) 

i 

= ^2 Yl < ' S ti S a,i+R +h - C -') ( 6 ) 

i 

G bb (R) = ^J2K S b ,i+R + h - c -) ( ? ) 

i 

Gab(R) = 4- h.c.) 

i 

= + M (8) 
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As (S+S~ +R + S~S+ +R ) = 2 (SfSf +R + S?S? +R ), they 
correspond to the spin correlations in the x — y piane. 

Adding the spin-spin correlations along the z axis 
(which are diagonal ąuantities) to the correlations in the 
xy piane that were obtained through Green functions, we 
obtain the complete spin-spin correlations. For example, 

S aa (R) = -^J2( S a,i-Sa,i+R). (9) 

i 

Similar definitions hołd for correlations S bb (R) between 
the b spins and for correlations S ab (R) between bosons 
and spins. We will denote by S tot = JF(S+ Sb.i) the 
total spin, or total magnetization, of the system, which 
is expressed as a sum of spin correlations functions 

‘Stot = E [S aa (R) + 2S ab (R) + S bb (R)} 

R 

III. PHASE DIAGRAM 

We first show the phase diagram in the (t/U,p/U) 
piane at T = 0 for a fixed value of V/U = 0.05 (Fig. 
0. For smali values of V/U, some of the transitions in 
the system are predicted to be first order 0 - At T = 0, 
using a canonical simulation, the Chemical potential is 
given by the energy difference p(N a ) = E{N a + 1) — 
E(N a ) where N a is the number of bosons. This allows to 
draw the boundaries of the nth Mott lobe with n bosons 
of type a per site by measuring the energy of the system 
with N a = nL 2 ,nL 2 + 1, and nL 2 — 1 particles. In 
the t/U = 0 limit, an analytical calculation yields 0 ] 
= -3V/4, = U- F/4, = 2 U- F/4, 

where /i n ^ n - t-i is the value at which the density changes 
from n to n+ 1. In Fig. [2j we exhibit the p = 1 and p = 2 
insulating Mott lobes. Outside of these Mott lobes, the 
system is superfluid and Bose condensed, as we will show 
below. The first Mott lobe is paramagnetic and the rest 
of the phase diagram has ferromagnetic correlations of 
the bosons. 

Compared to other studies of spin-1/2 bosons or 
mixtures of particles [l9}, the Mott lobes are larger. As 
expected, the presence of the Kondo interaction favours 
insulating behaviour: The tip of the p = 1 Mott lobe 
is located around t/U ~ 0.08 for V/U = 0.05 whereas 
it is located around t/U ~ 0.06 for F = 0. As V/U is 
increased up to 0.25, the tip shifts further to t/U ~ 0.10 
(see Fig. m- This robustness is not surprising sińce, for 
p = 1, the Mott gap is eąual to U + F/2 0 in the t = 0 
limit. 

Analysing the Green functions at p = 1, we observe 
that they all decrease rapidly to zero with distance in 
the Mott phase and there is no phase coherence (FigO 
(a)). This is expected for the one-body Green functions 
but is also the case for the anticorrelated Green functions. 
In the t = 0 limit, to minimize the AF energy between 
the spins and bosons, the magnetic moment of the boson 
forms a singlet with the spin located at the same site. As 



Figurę 2: Phase diagram for several system sizes L. The 
insulating phase at p = 1 is paramagnetic while the other 
phases are ferromagnetic. The dashed linę indicates the cut 
in Figurę [TO] 


these singlets are formed, there is a uniąue Mott State in 
the t = 0 limit and there is no possibility to exchange 
bosons of different spins 0] . This behaviour is probably 
maintained throughout the p = 1 Mott lobe, even at 
t / 0. This is verified by the value of G a b(0) ~ —0.5 
which shows the on-site antiferromagnetic correlation of 
the magnetic moments and will be confirmed below by 
direct measurements of the magnetic correlations. 

In the superfluid phase at p = 1 (Fig [3] (b)), on the 
other hand, all the Green functions show long rangę 
order. The long rangę order of the one-body Green 
function G ajCr shows that the system is Bose condensed. 
The non zero value of G aa shows that, in addition to the 
individual movement of the particles, exchange moves 
are important degrees of freedom for this phase. G bb 
is non zero, which shows that the spins are correlated. 
This correlation is mediated by the movement of the 
bosons as the spins are not directly linked to each other. 
This is confirmed by the observation that the boson- 
spin Green function, G ab , is also non zero. While 
G aa and G bb are positive, which signals ferromagnetic 
behaviour, G ab is negative which is expected sińce the 
coupling between bosons and spins is AF. The picture 
that emerges from these results is that the bosons and 
the spins form ferromagnetic phases but that these two 
species are coupled in an antiferromagnetic way (FigHJ). 

For p = 2, the Mott phase behaves differently from 
p = 1 (Fig. [5] (a)). The individual movement of 
particles are still suppressed: G a ,a goes to zero with 
distance. However exchange movements are present and, 
conseąuently, there are couplings between the spins as 
is seen from the anticorrelated Green functions taking 
finite values at large distances. This can be understood 
by noting that the ground State in the t = 0 limit is not 
uniąue [141 ]. This degeneracy will be lifted by a non zero 
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Figurę 3: The one-body and the anticorrelated Green 
functions versus distance R in the Mott insulator (a) and 
superfłuid phase (b) for p = 1. (a) all movement is suppressed 
as the particles form singlets. (b) movement of individual 
particles as well as anticorrelated exchanges are observed. 
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Figurę 4: The magnetic order that appears in the system. 
While the bosons and the spins form ferromagnetic layers, 
these two layers are coupled antiferromagnetically. 

hopping term and give a ground state with ferromagnetic 
correlations between the bosons. The magnetic order 
present in this phase is then similar to the one observed 
in the superfłuid phase. 

The p = 2 superfłuid phase (Fig. [5] (b)) shows the 
same qualitative behaviour as the p = 1 SF phase. 
However the dominant behaviour in this case is the 
anticorrelated movements of bosons whereas they were 
individual movements of bosons at p = 1. This is 
expected sińce, in a strongly correlated system, particles 
can move with a partner while individual movement is 
suppressed. Of course, as the density increases, the 
correlations become morę prominent. 


IV. QUANTUM PHASE TRANSITIONS 

We now analyse the naturę of the phase transitions 
between the Mott and the superfłuid phases by examining 
the behaviour of the superfłuid density. In the Mott 
lobe p s = 0; as t/U is increased at fixed p = 1 
and V/U = 0.05, we observe a seemingly continuous 
transition from the MI to the SF (Fig. [OJ. We fit the 
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Figurę 5: The one-body and the anticorrelated Green 
functions versus distance R in the Mott insulator (a) and 
superfłuid (b) phases for p = 2. (a) Individual movement are 
suppressed but anticorrelated moves are possible in the Mott 
phase. (b) All kind of movements are present in the superfłuid 
phase. 


curves near the transition with a form p s oc ( t/U — 
t/U c f (Fig. [6j and found /3 cs 0.33. It is however 
difficult to distinguish a continuous transition with a 
smali fi coefłicient (which gives a large slope close to 
the transition) from a discontinuous one on such finite 
size Systems. The condensate density p(k = 0) shows a 
similar behaviour (Fig. [TJ, as expected in two dimensions 
in the zero temperaturę limit where both p s and p{k = 0) 
are non zero in the superfłuid phase. A continuous 
transition is in disagreement with the MF analysis from 
[l4j that predicts a first order, discontinuous, transition 
when V/U < 0.1, which is the case here. To confirm 
this result, we calculated p{p) in the canonical ensemble 
and did not find a negative compressibility region which 
would have been a elear sign of a discontinuous transition 
[20( |. We also observe no sign of a discontinuity at the tip 
of the p = 1 Mott lobe, if there is one, as predicted in [1J] , 
it is too smali to discern for the system sizes accessible to 
us. In Fig. [HJ we show the evolution of this behaviour for 
a given size as t/U is increased. As V/U becomes larger 
the behaviour becomes smoother and the transition still 
appears continuous. 

In Fig. [9] we show p s as a function of t/U for the p = 2 
Mott-SF transition and do not observe a discontinuity 
between the two phases: the transition is continuous for 
all V/U. 

In figurę 1101 we examine the dependence of p and p s 
on the Chemical potential p along the dashed linę in 
Fig. [2] We observe the conventional incompressible Mott 
plateaux where dp/dp = 0. The superfłuid density goes 
to zero continuously as these plateaux are approached, 
showing that the transition for the p = 1 Mott is, as 
mentioned earlier, continuous. p(p) is also continuous 
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Figurę 6: Superfłuid density, p a , versus t/U in the first Mott 
lobe. We observe a continuous behavior for p s . The blue curve 
with triangles represents a linear extrapolation to L = oo of 
the results obtained for the three sizes L = 6, 8, 10. The 
curves for L — 10 and the extrapolated one have been fitted 
with a power law behavior p s oc {t/U — t/U c )^ (dashed lines). 
In both cases, we found t/U c — 0.79 and f3 ~ 0.33. 



Figurę 7: Condensate fraction, p(k = 0), versus t/U in 
the first Mott lobe. The shape is very similar to figurę [6] 
confirming the order of the transition. 

and has a positive slope, there is thus no sign of a phase 
separation close to the transition. 


V. MAGNETIC PROPERTIES 

We now study in morę detail the magnetic properties 
of the system. We plot S a p(L/ 2) at the largest possible 
distance L/2. For a sufficiently large system this 
converges to the sąuare of the magnetization of a given 
layer. We also study the on-site correlation to see if 



Figurę 8: Superfłuid density, p a , versus t/U in the first Mott 
lobe, for different ratio V/U. As we increase this ratio, p a 
becomes smoother, compared to Fig. [6] and the transition 
remains second order. 



Figurę 9: Superfłuid density, p„, versus t/U in the second 
Mott lobe. As for the p = 1 case (Fig. [6}, we do not observe 
any discontinuity in p a and the transition is second-order. 


singlets are formed between the a and b layers, and the 
total spin of the system S% ot /L A . 

In Fig. [Tl] we plot these ąuantities as functions of 
t/U for p = 1. We observe that in the Mott phase, the 
magnetic correlations are always zero. S a b{ 0) = —3/4 
signals the formation of a singlet. As the a and b 
particles form a singlet on each site, the absence of 
magnetic correlations between sites is reasonable. In the 
superfłuid phase, on the other hand, we observe that 
we no longer have a singlet phase as S a b{ 0) departs 
from the value observed in the Mott phase. We 
also observe, as anticipated earlier, that ferromagnetic 
correlations develop between the a bosons, between 
the b spins and antiferromagnetic correlations persist 
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Figurę 10: Cut of figurę [2] for t/U = 0.03, showing the total 
density, p, and the superfluid density, p s , as functions of p/U. 
Ali the transitions are continuous. 


between the two types of particles (corresponding to 
the positive values of S aa (L/ 2) and S b b(L/ 2) and the 
negative value of S a b(L/ 2)). Deep in the superfluid, 
the magnetic correlations take their maximum possible 
value \S a p(L/2)\ —> 1/4. It should be remarked that 
the correlation between the b spins is mediated by the 
itinerant a bosons, as there are no direct connections 
between the spins themselves. This is sinrilar to the 
coupling between localized spins provided by the RKKY 
interaction f22l [fki! in fermionie Systems, although it is 
always ferromagnetic in our case. Within error bars, 
we have S aa (L/ 2) = S bb (L/ 2) = -S ab {L/ 2) and, 
accordingly, the value of the total spin S/ ot /L 4 is zero. 
This was predicted in reference Ul in the high and Iow 
t/U lirnits but we see here that this seems to be the 
case also for intermediate values. We then have very 
different magnetic behaviour (independent singlets in the 
Mott, magnetic order in the SF) with the same Stot- 
One should remarks that, whereas the antiferromagnetic 
correlations between bosons and spins are imposed by the 
Hamiltonian and were expected, the ferromagnetism of 
the bosons layer appears spontaneously. There is no term 
that directly favors the development of FM correlations 
between bosons compared to other spin textures. 

Increasing V, we observe the same qualitative 
behaviour with some quantitative changes (Fig. 1121) . 
As mentioned earlier, the Mott-SF transition is shifted 
towards lower values of U as the Kondo interaction 
is added to the repulsion between particles, which is 
visible in the behaviour of p s . The appearance of the 
superfluidity and of the magnetic correlations is once 
again simultaneous and corresponds to the disappearance 
of the spin singlets. Finally we remark that the magnetic 
correlations tend to their maximum values but that those 
will be reached for much larger values of t/U. This 
is understandable as the singlets are morę diffleult to 



Figurę 11: Magnetic correlations versus t/U for p = 1 and 
V/U = 0.05. The dashed linę marks the transition between 
the Mott and the superfluid phases. There is no magnetic 
order in the Mott phase where the a and b spins form singlets, 
as shown by the values of S a b(0). The superfluid phase 
shows the magnetic behaviour depicted in Fig[4j intra-species 
ferromagnetic correlations and inter-species antiferromagnetic 
correlations. At large t, the magnetic correlations tend to 
their maximum values. 


break at large V which makes it morę diffleult to develop 
intersite correlations. 



Figurę 12: Magnetic correlations and p s versus t/U for p = 1 
and V/U = 0.25. The physical behaviour is the same as in the 
V/U = 0.05 case (Fig. 1111) but it is morę difhcult to establish 
the magnetic correlations with this larger value of V. 

For p = 2, deep in the Mott phase, the a spins 
located on the same site form a total s a s=a 1 moment, 
which gives 5 aa (0) = s a {s a + 1) ~ 2. This spin is 
then coupled antiferromagnetically to a b spin which is 
shown by the value of S ab ( 0) ~ — 1 , giving a total spin- 
1/2 and, consequently, two degenerate States on each 
site (Fig. fl3l) . The kinetic term lifts the degeneracy 
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between these States and we obtain the magnetic order 
which is observed in the SF regions and also even in 
the Mott lobe. Analytically 0 > it was predicted that 
S bb (R 0) = 1/36 ~ 0.0278. For the largest value of the 
interaction used in our simulations U = lOOt, we observe 
Sbb{L/ 2) = 0.028 ± 0.001 and, as can be seen in Fig. [T31 
we reached the regime where Sbb{L/ 2) saturates at this 
non zero value. 

In the superfluid phase, the system behaves very much 
as for p = 1. 5 oo (0) and S a b( 0) depart from their Mott 
phase value and increase slightly. S aa (L/ 2), S/ & (L/2), 
and S a b(L/ 2) go to their extreme possible values 1, 1/4 
and -1/2, respectively. 

It is predicted 0 that 8/ 2 ot = (L 2 /2)(L 2 /2+l) ss L 4 /4 
in the strong and weak coupling limits. We observe that 
5/ ot /L 4 always takes a value compatible with 1/4, for 
any value of t/U. As for p = 1, we observe two different 
behaviours for the same common value of Slt/L 4 . 



Figurę 13: Magnetic correlations versus t/U for p = 2 and 
V/U = 0.05. The dashed linę marks the transition between 
the Mott and the superfluid phases. The same magnetic order 
is present in the Mott and the superfluid phases although 
there are two different limiting regimes with stronger (for 
t/U 3> 1) or weaker (for t/U <§; 1) b spin correlations. The 
magnetic behaviour is the one sketched in Fig. [4] The total 
spin seems constant Sfot/L 4 ~ 1/4. 


VI. THERMAL EFFECTS 

We analysed the behaviour of the observed phases at 
finite temperaturę. First we looked at the superfluid 
density to determine the extent of the superfluid phase 
as the temperaturę is increased. The thermal phase 
transition between the superfluid phase and the normal 
liąuid is of the BKT type. We performed different 
finite size analyses to determine the critical temperaturę 
T c at which p a becomes zero. First we used linear 
extrapolations of p s as a function of 1/L for different 


values of 1/T. Then we used Nelson and Kosterlitz’s 
result 0 p s {T c ) = kT/nt to calculate the temperaturę 
T C (L) at which our curves intersect kT/nT before looking 
at the 1/L —>• 0 extrapolation (Fig. [Ul) . Finally we used 
the recently proposed method by Hsieh et al. 0. Ali 
three methods gave similar results for the system sizes 
accessible to us. 



Figurę 14: Superfluid density, p s , versus T for several sizes L 
with t/U = 1 and V/U = 0.05. We calculate the intersection 
between the curves and kT/nt to obtain T C (L). We then 
extrapolate linearly (see inset) the value of T C (L) to the large 
size limit to obtain the estimate of the transition temperaturę 
shown in Fig. 1161 

There is no transition between the Mott phase and 
the normal liąuid, as the Mott phase exists only at 
zero temperaturę, strictly speaking. We calculated the 
fluctuation of the number of particles k = (n 2 /) — ( n a ,i ) 2 
which exhibits a plateau at smali kT before increasing 
at higher kT. We identify the crossover temperaturę 
between the Mott and the liąuid behaviour as the 
temperaturę where k departs from this Iow T value by 
morę than 5%. We checked that this definition is valid 
by comparing with a measure of p(p) and finding the T 
at which the Mott “plateaux” disappear in p(p). 

Putting these results together, we obtain the phase 
diagrams shown in Fig.[l6]for p = 1 and p = 2. We placed 
at T = 0 the point of the ąuantum phase transition 
observed in Sec. ED 

We now turn to the magnetic behaviour at finite 
temperaturę. In the Mott phase at p = 1, there is no 
magnetic order at T = 0 and this behaviour persists 
at finite temperaturę. In the Mott phase at p = 2, 
the magnetic couplings that lead to a ferromagnetic 
phase are weak and the magnetic order disappears for 
Iow temperatures kT ~ 0.3 on an L = 6 system. As 
the system is well described, in the Mott phase, by a 
Heisenberg model, we do not expect to observe magnetic 
order at finite temperaturę in the thermodynamic limit. 

The behaviour in the superfluid regime is much morę 



































Figurę 15: The local density fłuctuation k = (n\ — (n a ^) 2 
versus kT for different temperatures. We see that k is 
alrnost constant at Iow kT in the Mott region. We defined 
the crossover temperaturę between the Mott and the liąuid 
behaviour as the temperaturę for which k departs from its 
Iow T value by morę than 5%. 



Figurę 16: The phase diagrams for p = 1 and p = 2 as a 
function of t/U and kT/U at V/U = 0.05. There are three 
different regions in each diagram: a Mott region (MI), the 
normal Bose liąuid (NBL) and the superfłuid region (SF). 
The SF-NBL limit is a BKT transition. The MI-NBL limit is 
a crossover between weakly and highly compressible regimes. 


interesting. What we observe is that the magnetic 
ordering is reinforced when we increase the temperaturę 
from T = 0. We hnd that the total magnetisation Ąot 
increases before decreasing again and reaching zero in the 
high temperaturę regime (see Fig. [171) . The increase of 
the total magnetisation is due to a strong decrease of the 
boson-spin S a b correlations. This is easily understood 
as the coupling V between bosons and spins takes a 
smali value (V/U = 0.05). Thermal excitations break 


Figurę 17: Magnetisation of the whole system (S to t) and 
sums over all distances of boson-boson (S aa ), spin-spin (Sbb) 
and boson-spin ( S a b ) magnetic correlations for p = 1 in the 
superfłuid phase. The total magnetisation Stot is zero at 
T = 0 but rises when T is increased. As the bosons decouple 
from the spin, when kT > V, and no longer form singlets with 
those, they develop ferromagnetic correlations. Top panel: 
p — 1, bottom panel: p = 2. 



Figurę 18: Extrapolation of the magnetization as a function of 
the size. At high temperaturę (kT = 4.55) the magnetization 
extrapolates to zero. At the temperatures where we observed 
the peaks of the moment in Fig. [TT] the magnetization does 
not seem to extrapolate to zero. The data for p = 2, kT = 0.9 
have been divided by 2 for better visibility. 

the correlations between spins and bosons. The AF 
correlation between those two species disappears and the 
spins are then disordered as the bosons no longer mediate 
an intersite coupling. As the bosons become independent 
of the spins, they form a FM superfłuid with a larger 
total magnetisation. Once again, the mechanism for this 
FM ordering is not elear but it is obviously mediated 
by the hopping of bosons, the only effect that can couple 
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distant particles in this system. As the hopping t is larger 
than V these remaining FM correlations disappear only 
at larger temperatures. For example, for p = 1 (Fig. [TO 
the AF correlations S a b and the spin-spin correlations 
Sbb have almost disappeared for kT = 1 whereas the FM 
correlations S aa become negligible only for kT > 3. 

A non zero magnetisation should not be present at 
finite temperaturę in two dimensions as it contradicts 
the Mermin-Wagner theorem [23]. We looked at the 
evolution of the magnetisation with sizes for different 
temperatures (Fig. [TO]) to check if it decays to zero. 

At high temperaturę (kT = 4.55) the total 

magnetisation indeed goes to zero. With no correlations 
between sites, S^ ot scales as L 1 2 * , as the only remaining 
contributions are on site terms. Hence the observed 
behaviour where S^ ot /L 4 oc 1/L 2 . However, at 
intermediate temperaturę, close to the maxima of iS 2 ot 
observed in Fig. [TT] we do not find a elear decay of 
the magnetisation with size. Data obtained for L = 12 
and L = 14 have large error bars and it is difficult to 
draw a conclusion for the thermodynamic limit, but the 
behaviour does not seem to correspond to an exponential 
decay of the magnetic correlations with distance. 

Our interpretation of these data is that, in the 
superfłuid phase, we have a quasi-long rangę order of the 
different Green functions. We are then also expecting 
a quasi-long rangę order for the magnetic correlations 
between the bosons as they are directly related to the 
anticorrelated Green functions. This would explain the 
surprisingly large values of S 2 ot on our smali size Systems. 
A similar behaviour was found in another spin 1/2 bosons 
model at finite temperaturę [28j . 

VII. CONCLUSION 

We have studied a Bosonic Kondo-Hubbard model 
with an AF interaction between spin-1/2 bosons and 
hxed spins. We have drawn the phase diagram of 
the system and found that the presence of the Kondo 
interaction with the spins facilitates the localisation of 
the particles into Mott phases. We have shown that 
exchange moves are taking place in Mott phases with 
density larger than one. Studying the naturę of the 
phase transition, we have always observed continuous 
transitions, contrary to the MF prediction [l4| that a 
discontinuous transition is present at the tip of the p = 1 
Mott lobe for Iow enough V. 

The magnetic properties of the system are particularly 
interesting. At zero temperaturę, the total magnetization 


of the system is always constant but different behaviours 
can nevertheless be observed. In the p = 1 Mott phase, 
we have observed on-site singlets between the bosons and 
the spins, with no long rangę order. On the contrary, in 
the superfłuid, we have a FM order of the bosons and the 
spins and an AF order between them. At p = 2 we always 
have this same magnetic behaviour for all interactions 
but with two limiting regimes in the Mott and superfłuid 
phases. Notably, we observed the very smali value of 
the spin-spin correlations at large U that was predicted 
analytically. 

At finite temperaturę, we determined the boundary of 
the superfłuid phase and found the crossover temperaturę 
between the Mott and liquid regions. Morę interestingly, 
we found that, in the superfłuid phase, the total 
magnetisation is inereased due to the fact that the bosons 
decouple from the spins. This is unexpected and is 
certainly not present in the thermodynamic limit but 
should be observed on finite size Systems. 

The underlying physics of the boson Kondo 
Hamiltonian studied here has significant similarities, 
but also several differences, from the fermionie case. For 
fermions at commensurate hlling (p = 1), although there 
is a singlet-antiferromagnetic phase transition, both 
magnetic phases are insulating, whereas bosons at weak 
coupling are superfłuid. The naturę of the magnetic 
order is also somewhat different. In the fermionie case, 
ordering of the local spins separated by a distance r 
is mediated by a Rudermann-Kittel-Kasuya-Yoshida 
interaction which has a modulation cos(k F • r) where 
k F = ( 7r, 7r) is the Fermi wavevector at p = 1. In 
contrast, the order in the bosonic case studied here is 
ferromagnetic. There is much past and current interest 
in the Kondo Hamiltonian for fermions for various sorts 
of dilution and randomness in order to model novel 
quantum phase transitions and also Chemical doping of 
heavy fermion materials M. It would be interesting to 
study analogous effects in the boson-Kondo Hamiltonian. 
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